Flow rule of dense granular flows down a rough incline 
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We present experimental findings on the flow rule for granular flows on a rough inclined plane 
using various materials including sand and glass beads of various sizes and four types of copper 
particles with different shapes. We characterize the materials by measuring h s (the thickness at 
which the flow subsides) as a function of the plane inclination 9 on various surfaces. Measuring the 
surface velocity u of the flow as a function of flow thickness h, we find that for sand and glass beads 
the Pouliquen flow rule u/^/gh ~ /3 h/h s provides reasonable but not perfect collapse of the u(h) 
curves measured for various 9 and mean particle diameter d. Improved collapse is obtained for sand 
and glass beads by using a recently proposed scaling of the form u/^fgh = f3 ■ h tan 2 9/h s tan 2 9\ 
where 9\ is the angle at which the h 3 (9) curves diverge. Measuring the slope ft for ten different 
sizes of sand and glass beads, we find a systematic, strong increase of /3 with the divergence angle 
#i of h 3 . The copper materials with different shapes are not well described by either flow rule with 
u ~ h 3 ' 2 . 

PACS numbers: 47.57.Gc, 45.70.-n 



I. INTRODUCTION 

Granular flow on a rough inclined plane is an important 
system with which to learn about the basic rules of the 
dynamics of granular materials p], 0, 0, 0, 0, 0, 0, 0, 0] . 
Despite intensive study, the fundamental features of such 
flows are still incompletely understood (for reviews see: 
[UlEl EH)- The majority of laboratory experiments re- 
port on the flow properties in narrow channels (quasi-2D 
geometry) where the velocity can be measured as a func- 
tion of depth by directly vie wing gr ain motion through 
the sidewalls @, @, 0, i, H, M, GJ, E3, d ■ In this config- 
uration, however, the effect of friction with the confining 
vertical walls is important [HI, HH, E3] > an d remains 
a determining force for thick flows (flow on a pile) even 
in wider channels [lj], EH ■ 

For thin flows in wide channels measuring the depth 
dependence of the flow velocity (far from the side walls) 
is far more difficult. To characterize the basic features of 
granular flows in this configuration, the surface velocity 
u or the depth averaged velocity U can be measured as a 
function of the flow thickness h. The depth averaged flow 
velocity U, inferred from the front velocity of the gran- 
ular layer, was systematically measured by Pouliquen as 
a function of the flow thickness h for glass beads over a 
range of plane inclinations 9 Q. The U(h) curves mea- 
sured at different values of 9 collapsed when using the 
scaling law U/^/gh = (3 h/h s — 7 (where h s corresponds 
to the thickness where the flow subsides) giving rise to a 
general flow rule, denoted the "Pouliquen flow rule," for 
glass beads with various sizes and for which f3 ft 0.14 and 
7 ft 0. It was subsequently reported [|[ that the same 
scaling collapsed the U (h) curves for sand with one par- 
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ticular size of d = 0.8 mm. The slope for the sand data, 
(3 ft 0.65, was considerably larger than for glass beads 
and 7 ft 0.77. This quantitative difference in the flow 
rule was used to explain complex dynamical phenomena, 
such as waves Q or avalanche propagation [l9| ■ 

It is of considerable interest to determine the robust- 
ness of the Pouliquen flow rule (PFR) for different flow 
conditions including particle diameter, relative surface 
roughness, and particle shape. A further consideration is 
whether the flow rule is sensitive to measuring the front 
velocity as compared to measuring the surface velocity. 
The former yields a better depth averaged velocity but 
is subject to saltating grains for faster flows which lim- 
ited the accuracy of the measurement to about 10 % 
and would not be applicable for gener al g ranular mate- 
rials subject to a fingering instability [20|. The surface 
velocity measurement is characteristic of the steady flow 
even for general granular media and avoids the accuracy 
limitations imposed by saltation, but can only be related 
to a depth-average by some assumption of the vertical 
velocity profile. Neither approach is ideal, being rather 
complementary as opposed to one being a priori better 
than the other. Our work establishes the utility and ro- 
bustness of using the surface velocity to determine the 
flow rule. If a flow rule is to be a useful measure of the 
state of granular flow on an incline, it should not be par- 
ticularly sensitive to the details of its determination. 

A recent theory by Jenkins plT | suggests a phenomeno- 
logical modification of the hydrodynamic equations for 
dense flows. According to the theory enduring contacts 
between grains forced by the shearing reduce the colli- 
sional rate of dissipation while continuing to transmit 
force and momentum. This assumption has several con- 
sequences, one of which is a modification of the Pouliquen 
scaling law by the inclusion of a tan 2 9 correction to 
the h/h s term. Rcplotting the Pouliquen data, Jenk- 
ins found a better collapse of the data using his modified 
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form, denoted here as the "Pouliqucn- Jenkins" flow rule 
(PJFR). The improvement of the collapse, however, was 
not definitive owing to the scatter in the velocity data 
and in the associated determination of the h s (9) . 

One of the purposes of a flow rule is to have a com- 
pact description of easily measurable quantities that rep- 
resents the subtle balances of stress and strain rate in 
a granular material, i.e., the granular rheology. Al- 
though the vertical velocity profile in a flowing granular 
layer has not been obtained experimentally, let alone the 
experimental determination of local stresses and strain 
rates, a general discussion about possible flow rhcolo- 
gies helps set a background for presenting empirical flow 
rules determined from experiment. In particular, the 
scaling of the velocity with layer thickness can be un- 
derstood by a consideration of bulk Bagnold rheology 
[EH2]- I n the theory of Bagnold, the shear stress varies 
with the shear rate 7 like <r xz ~ j 2 . This relationship 
is based on the following assumptions. The transport 
of the x component of momentum in the z direction oc- 
curs through collisions whose rate depends on the ve- 
locity gradient 7. Similarly the momentum transfer per 
collision scales linearly with 7 leading to the quadratic 
dependence between stress and strain. With a linear de- 
pendence of shear stress on the vertical coordinate z, this 
leads to a vertical variation of the down-plane velocity of 



u(z) ~ h 3 ' 2 
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Thus, the surface ve- 



{(h~z)/hf 2 

locity u = u(h) ~ b?l 2 so that the scaling u/y/gh versus 
h (suitably corrected for inclination angle) should yield 
straight lines with zero intercept. Such scaling was re- 
ported for experiments 0, [23| using glass spheres and for 
numerical simulations of idealized spherical particles [22| . 
Also, deviations from this law towards a linear velocity 
profile were reported in experimental [23| and numerical 
[2J] studies for thin flows. For a particular flow profile, 
the surface velocity u and the depth averaged velocity 
U are related by a constant factor. Thus, there is no a 
priori reason to prefer one over the other. Although the 
interior velocity profile far from sidewalls has not been 
measured to our knowledge, the scaling u ~ h 3 / 2 is indi- 
rect support for the Bagnold flow rheology. The degree 
to which such scaling fails, therefore, would appear to 
call for modification of the assumptions leading to the 
Bagnold rheology. We will sec in this paper how well the 
Bagnold-based rheology applies to a range of different 
granular materials. 

In the present work we investigate the flow prop- 
erties of 14 different materials by measuring the sur- 
face velocity u as a function of flow thickness h. Be- 
cause of our measurement methods, the statistical un- 
certainty in our data is considerably less than in previ- 
ous studies 0, allowing for a more detailed and quan- 
titative evaluation of different flow rule searings. We 
find that scaling the surface flow velocity by \fgh and 
the flow thickness by h s assuming the Pouliqucn flow 
rule provides reasonable but not perfectly accurate col- 
lapse of the u{h) curves taken at various plane inclina- 



tions measured in a wider range of the main control pa- 
rameters of grain size, plane inclination angle, surface 
roughness and flow thickness compared to earlier studies 
[3, HI [HI- Improved collapse is obtained for sand and 
glass beads using the modified Pouliquen-Jenkins scal- 
ing law u/y/gE = (i ■ h tan 2 6/h s tan 2 9\ where the factor 
tan 2 6 is supported by a recent theory [2l[ . For glass 
beads the straight lines of the scaled curves support the 
Bagnold rheology. For sand, although the data are well 
collapsed by the scaling, the curves arc slightly concave 
downward suggesting high-order corrections in h beyond 
the simple Bagnold result. We show that the slope (3 of 
the master curve for the sand/glass-bead materials (ob- 
tained for each material) strongly increases with tan#i 
or tan(9 r where 9\ and 9 r are the angles where h s (9) di- 
verges and the bulk angle of repose, respectively. The 
similarities and differences of our experimental approach 
compared to other experimental measurements of flow 
rules @, are discussed in detail. 

In contrast to the relatively simple and understandable 
data obtained for glass beads and for sand, the behavior 
of flowing copper particles is more complex and a sim- 
ple Bagnold interpretation works quite poorly in describ- 
ing the relationship between surface velocity and layer 
height. Indeed, the the scaling of u with h is closer to 
u ~ h 1 / 2 than to the Bagnold form u ~ h?l 2 . Never- 
theless, the angle correction using h s (8) (Pouliquen flow 
rule) or h s (6) / tan 2 6 (Pouliquen-Jenkins flow rule) ap- 
pears to work pretty well with the latter again providing 
better overall data collapse. 



II. EXPERIMENT 

The experimental measurements presented in this pa- 
per were performed in two different setups. The first 
apparatus was described in detail elsewhere [25| and con- 
sisted of a glass plate with dimensions 230 cm x 15 cm 
(see Fig. [TJ . The leftmost 40 cm of the plate served as 



hopper 




FIG. 1: (Color online) Schematic of the experimental setup 
used to measure the surface velocity and height. The whole 
system could be rotated (together with laser and camera) to 
set an arbitrary inclination angle 9. 
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the bottom of the hopper. The surface of the remaining 
part (190 cm) of the glass plate was typically covered 
with sandpaper that was glued to the surface and had a 
roughness of R = 0.19 mm (grit 80) which provided an 
extremely durable uniform roughness to the plate. Other 
values of plate roughness were studied using different grit 
sandpaper (and a few measurements with 0.4 mm sand 
glued to the glass plate) to explore the systematic de- 
pendence of our results on relative roughness compared 
to grain size. The plate together with the hopper could 
be tilted, enabling us to set an arbitrary inclination angle 
6. The flow was characterized by measuring the surface 
velocity u as a function of the thickness h in the station- 
ary dense flow regime at a location x = 155 cm below 
the hopper gate, sufficiently far downstream to have es- 
tablished a steady state [25|]. This first apparatus could 
be tilted back and forth to recharge the hopper, facil- 
itating the accumulation of the large amounts of data 
reported here. Because the system was closed in a cylin- 
drical tube, precise measurements of h s {9) were difficult 
and were performed in a second apparatus. 

The second apparatus used to measure h s {9) consisted 
of a wider plate having dimensions 227 cm x 40 cm and 
was covered with the same sandpaper used to cover the 
surface of the flow channel. The system was not con- 
fined from the top. This wider, non-enclosed channel 
allowed for very precise measurements of the layer height 
in a rapid manner. The procedure was to throw the 
grains onto the plane and allow a uniform layer thick- 
ness to form by letting the flow subside. The granular 
material was swept from a 2 m long area and its volume 
was measured accurately, yielding a very precise and re- 
peatable measurement of the mean layer thickness h s (8). 
This method averages out the spatial variations in h s , 
the amplitude of which was also estimated by measuring 
the displacement of a projected laser sheet. At lower 8, 
the height variations were typically less than ±5% of h s , 
but became larger at higher plane inclinations where h s 
became less than 5c?. Because the majority of the data 
on the flow properties were measured at plane inclina- 
tions corresponding to these relatively lower values of h a , 
it is important to get an accurate measure of h s . The 
repeatability of the measurement also depended on the 
plane inclination, but in this case relative variations de- 
creased with increasing 8. The data points fell onto the 
same curve within an error of ±5% for tan 8 / tan 8 r > 1.1. 
When approaching 8 r the measurements became less ac- 
curate with the rapid increase of h s leading to a ±12% 
uncertainty of the data points for tan 8/ tan6> r < 1.1. 

Uncertainties arising from slightly nonuniform thick- 
ness near the walls were also estimated. We observed a 
boundary layer W where the layer thickness was slightly 
larger than elsewhere. The width of the boundary layer 
was W < 1 cm for tan 9/ tan r > 1.1 and somewhat 
larger W < 5 cm for smaller 8. The effect of the bound- 
ary layer results in a slight overcstimation of h s corre- 
sponding to about 2% for tan 8/ tan6> r > 1.1 and 5% for 
smaller 8. To reduce the effect of the boundary layer we 



removed part of the excess material near the boundary, 
and we estimate that the finally measured value of h s is 
overestimated by less than 1% owing to the effect of the 
lateral boundaries. 

A possible concern regarding using one apparatus to 
measure h s {8) and another to measure u and h in the 
flowing state is that the lateral boundary effects might 
be different, leading to possible discrepancies in the mea- 
surements. To that end we measured the flow profile in 
the narrow channel, as presented below, and found that 
the flow was uniform over the central 80% of the nar- 
row channel. Because our measurements of u and h were 
taken in the center of the narrow channel, we conclude 
that no significant differences arise from using different 
channels for the static and dynamics measurements, re- 
spectively. Further, because of the limitations of each 
system, the amount of data we obtained would not have 
been feasible using one or the other of our experimental 
setups. 

Four types of granular materials were used. The first 
set consisted of sand particles from the same origin but 
sorted into four different sizes. For example, the finest 
sample was obtained by sifting the sand with 100 and 
300 /im sieves. We designate this distribution as having 
a mean of d = 0.2 mm and a standard deviation of 0.05 
mm. According to this notation the four sets of sand 
correspond to sizes d = 0.2 ± 0.05 mm, d = 0.4 ± 0.05 
mm, d = 0.6 ± 0.05 mm and d = 0.85 ± 0.08 mm while 
the mean particle density was p sa nd = 2.6 g/cm 3 . The 
fifth sample of sand originated from the Kelso dunes and 
was well sorted with a size distribution of d = 0.2 ± 0.05 
mm. This sand is peculiar in that it emits sound when 
sheared. The Kelso dune is known to be an example of 
"booming sand dunes" [HJ. We also used commercial 




FIG. 2: (Color online) Microscopic images of the copper 
particles with d — 160 ± 50 /im and with packing fractions r\: 
a) 0.25, b) 0.33, c) 0.5 and d) 0.63. 

glass beads (Cathapote) with sizes d = 0.18 ± 0.05 mm, 
d = 0.36 ± 0.05 mm, d = 0.51 ± 0.05 mm and d = 0.72 ± 
0.08 mm and mean particle density of p g i aS s = 2.4 g/cm 3 . 
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One sample of the d — 0.51 ± 0.05 mm glass beads was 
carefully washed. For that sample, we observed a slight 
change in the flow properties as well as in the value of h s 
compared to an unwashed sample with the same d and, 
thus, we report these data as an additional case. The last 
type of material consisted of copper particles with a mean 
size of d = 0.16±0.03 mm but with different shapes. The 
shape anisotropy of the four different samples of copper 
particles is characterized by the volume fraction r] (the 
ratio of the volume occupied by the particles and the 
total volume) of the material at rest, with values 0.63, 
0.5, 0.33, 0.25 and particle densities 8.7, 8.2, 7.6, 7.1 
g/cm 3 , respectively. The variation in r\ represents the 
strong change in the shape from spherical particles to 
very dendritic shapes with decreasing particle densities 
for the more dendritic shapes as well. Images of copper 
particles are shown in Figs. [2ja)-[2jd) where the strong 
variation in particle shapes is clearly seen. 

The surface flow velocity u was determined by analyz- 
ing high speed (8000 frames per second) video recordings. 
Individual particles make streaks in a space-time plot 
of intensity along one line of camera pixels aligned with 
the mean flow direction. An example of such a space- 
time image is shown in Fig. [2Kb) where the length of the 
line in the camera is L = 3.68 cm and the total time is 
T = 0.080 s. The streaks are generally oriented at some 




b. 


/<x\ 







FIG. 3: a) Space time plot showing particle streaks along a 
line oriented with the flow direction for 6 = 36.1° and H — 2 
cm. Dimensions of the image are 3.68 cm and 0.08 s. b) 
Two dimensional FFT in frequency /wave- number space of the 
image in a) which gives an accurate measure of the mean 
flow velocity as indicated by the solid line. The angle a is 
designated in each image. 



angle a in the image. Performing a fast Fourier trans- 
form (FFT) produces a line perpendicular to the streaks, 
see Fig. [31 which gives a measure of the mean surface ve- 
locity u = (L/T) tana. The thickness h of the flow was 
monitored by the translation of a laser spot that was pro- 
jected onto the surface of the plane at an angle of 4> = 20° 
in the xz plane (see Fig. [TJ. Other details regarding the 
measurement techniques can be found in [25| . 



III. RESULTS AND DISCUSSION 

The two measurements that determine the flow rule are 
the height of the layer when the flow stops h s {ff) and the 
dependence of the surface velocity on the layer height h. 
We first consider h s {9) for glass beads, sand and copper. 
We then present measurements of u as a function of h for 
sand and glass beads and the application of the flow rules 
of PFR and P JFR. Finally, we consider velocity data and 
flow rules for the copper material. 



A. Determination of h 3 (0) 

As seen in Figs.[3ta)-[lKc), h s increases rapidly with de- 
creasing 9 and diverges at 9\. The solid lines in Fig. 2] are 
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FIG. 4: (Color online) The thickness h s at which the flow 
subsides normalized by the grain diameter d as a function of 
the plane inclination angle 9 for a) sand and b) glass beads 
of various sizes and c) copper particles of various shapes (as 
indicated by the static volume fraction if). The grain diameter 
d is indicated and w designates the case of washed glass beads. 
The continuous lines are best fits according to the formula 
h s /d — B (tan 62— tan 6) /(tan 6 — tan#i). The resulting values 
of #1 are indicated for each material in Table [I] 

best fits to the formula h s /d = ^4(tan 6*2 — tan 9) / (tan 9 — 
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tan#i). a simple function which diverges at 9\ and goes 
to zero at 9 = 9 2 0, H, [HI H3] ■ The resulting values for 
the fitting parameters A, 9\, and 9 2 are indicated for each 
material in Table [J The bulk angle of repose 9 r was also 
measured for several materials by measuring the dynam- 
ics of a three dimensional sandpile under constant flux 
conditions. As material was added at a very small but 
uniform rate to the top of the pile, avalanches formed 
and propagated downward intermittently. The distribu- 
tion of the angle, observed directly after the avalanche 
stopped, was measured for hundreds of avalanches. The 
mean of this distribution was taken to be 9 r , the bulk 
angle of repose. The value of 9 r is very close to 9\ as 
indicated in Table [U 

The h s (9) curves are very similar for all four sand sam- 
ples originating from the same source [see Fig.0Ja)]. The 
fifth curve corresponding to the Kelso sand showed devi- 
ations from the other data at lower values of 9, yielding 
a somewhat smaller value for 9\. This difference is at- 
tributable to the more rounded shape for the Kelso sand 
as revealed in microscope images. The h s {9) curves for 
glass beads, however, formed two groups. Microscope 
images revealed that the two samples with smaller d con- 
tained a larger amount of non-spherical particles than the 
two sets with larger d. 
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TABLE I: The values of 8 r , 0i, di and A for sand and glass 
beads of size d and copper particles with d — 0.16 mm and 
volume fractions r\. K stands for the Kelso sand, and w de- 
notes the washed glass beads. 

The difference in shape may explain the slightly larger 
values of 9 1 and h s /d for the two samples with smaller d. 
The case of the 0.51 mm glass beads is also interesting 
in that washing the material with tap water resulted in a 
slightly smaller value of h s , which implies slightly reduced 
friction either with respect to the rough surface or be- 
tween individual grains. The reduction could have been 
caused by the elimination of non-spherical dust particles 
owing to washing the sample. The four samples of copper 



are nice examples of the effect of particle shape. A sys- 
tematic increase of 9\ and 9 r detected by changing shape 
anisotropy in the order of spherical beads (rj = 0.63), 
particles with irregular but rounded shapes (77 = 0.5) and 
the two sets of particles with very anisotropic dendritic 
shapes (rj = 0.25 and rj = 0.33). 

The influence of the boundary conditions can have a 
profound effect on the conditions of the granular flow. 
The usual no-slip boundary condition appropriate for a 
fluid is probably never completely satisfied for a granular 
flow and certainly depends on surface roughness. Fur- 
ther, the role of the surface in damping energy is only 
recently beginning to attract attention [281 ] and has not 
been considered in the context of granular flows on an 
incline. Thus, it is important to evaluate the dependence 
of our results on surface roughness and, in principle, on 
surface restitution coefficient. Although we do not con- 
sider the latter here, the systematic of a flow rule com- 
parison may depend on the damping properties of the 
surface which may help explain differences between flow 
on a soft felt surface, on a glass plate with glued on hard 
particles, or a hard surface covered with sandpaper. 

To study the dependence of our results on surface 
roughness, we measured the dependence of h s /d on plane 
roughness R, shown in Figs.[5Ja) and[S{b) for sand with 
d = 0.4 mm and glass beads with d = 0.36 mm on four 
different sandpapers with nominal roughness of R = 0.12 
mm, R = 0.19 mm, R = 0.43 mm and R = 0.69 mm 
(grits 120, 80, 40 and 24 respectively). For sand h s /d 
was also determined on a surface prepared by gluing one 
layer of the same grains onto the plate. For both sand 
and glass beads, Figs.^a) and[5]^b), a slight increase of 
h s /d is observed with increasing plane roughness. For 
the case of sand with d = 0.4 mm the curve measured 
on the surface prepared by gluing the same grains was 
the most similar to the curves taken on sandpaper with 
R = 0.12 mm or R = 0.19 mm, i.e., the surface friction 
for sandpaper is somewhat larger than when the surface 
is covered with sand glued to the surface. 

We determine the relative effect of surface roughness 
on the determination of h s by comparing data for sand 
and glass beads for different values of d and R. The value 
of h s /d increases as a function of R/d as shown for three 
values of tan0/tan0 r in Fig.OJc). At plane inclinations 
close to the bulk angle of repose 9 r , the curve seems to 
saturate [see the curve taken at tan 6/ tan 9 r = 1.1 in 
Fig. Efc)] but for larger plane inclinations, i.e., for thin- 
ner layers, a slight increase of h s /d is observed over the 
measured range of R/d. The increasing tendency of h s /d 
indicates that the effective friction near a rough surface 
increases slightly with increasing plane roughness. Near 
the rigid surface the particles have less freedom to rear- 
range so that in order to shear the medium has to dilate 
more [29| yielding a larger effective friction, compared 
to the case of the bulk material. The growing value of 
hs/d matches the overall tendency of the data reported in 
[lfj using monodisperse glass beads on surfaces prepared 
by gluing one layer of glass beads on a plate. We did 
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FIG. 5: (Color online) h 3 vs. 6 for a) sand with d — 0.4 
mm and b) glass beads with d = 0.36 mm for various values 
of surface roughness: surface covered by sandpaper with R — 
0.12 mm (x), R = 0.19 mm (*), R = 0.43 mm (o) and R = 
0.69 (o), surface covered by one layer of d — 0.4 mm sand 
particles (A), c) h s /d as a function of R/d for tan#/ tan# r = 
1.1 (x); 1.25 (*); and 1.4 (o) for sand and glass beads. 



not, however, find any significant height maximum cor- 
responding to a particular plane roughness reported in 
[l(J[3(|. Note that a stronger difference in h s /d was de- 
tected when the values measured on a solid rough surface 
(similar to our case) and on velvet cloth were compared 
[nil ]. As discussed above, this may be more a result of 
surface damping than surface roughness. 

The homogeneous dense- flow regime existed for moder- 
ate plane inclinations where tan 9/ tan 9\ was in the range 
1.1 — 1.45. According to our measurements (25[, the den- 
sity of the flow in this regime was decreased slightly with 
increasing 9 but was always larger than 0.8p s where p a 
is the close packed static density of the material, in ac- 
cordance with other experimental data 0, H| and with 
numerical simulations [22, n24, Kllf . 



B. Flow rule for glass beads and sand 

We next present measurements of flow velocity ob- 
tained using the space-time technique described above. 
To demonstrate that sidewall boundaries do not affect 
the velocity near the channel center, we consider the 
transverse velocity profiles shown in Fig. [5] for a hopper 
opening of H = 2 cm for several values of 9. The data 
show that friction with the smooth sidewalls is much less 




y (cm) 

FIG. 6: (Color online) Transverse velocity profiles of the flow 
for sand with d = 0.4 mm at a hopper opening of H = 2 cm. 

important than friction with the rough bottom plate so 
that the sidewalls only produce a lateral boundary layer 
at the edge of the channel with a characteristic thickness 
of 2-3 cm. Over the remaining 80% of the channel width, 
u is very constant. For determination of the flow rule, the 
velocity u and thickness h were measured at the channel 
center. 

We now consider u as a function of h for sand and glass 
beads, presented in a variety of forms to test both PFR 
and PJFR [p. In Figs. E^a) andE[b), we show u as 
a function of the flow thickness h for sand with d = 0.4 
mm and glass beads with d = 0.36 mm. In Figs. [Tfc) 
and [7Jd) , the same data are presented in dimcnsionlcss 
form accordingto the flow rule u/y/gE ~ h/h s suggested 
by Pouliquen [2|. For comparison, we include the curves 
measured by Pouliquen for glass beads with d = 0.5 mm 
and sand with d = 0.8 mm correcting for the dif- 
ference between depth averaged velocity U and surface 
velocity u that assumes a Bagnold velocity profile for 
which u = 1.67U. 

Our data cover a wider range of u and h than previ- 
ously measured, partly because of the smaller grain size, 
but also owing to the measurement technique. Namely, 
measuring the surface velocity in the stationary regime 
was much more straightforward for us than detecting the 
velocity of the front and thereby determining U . For the 
detection of the front velocity the difficulty was that in 
contrast to the simple monotonic increase of the height at 
the flow front (reported in 0]), in some cases and partic- 
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FIG. 8: (Color online) Dimensionless flow velocity u/ \fgh vs 
h/h s for sand and glass beads as a test of PFR. The grain 
diameter d is indicated and w designates the case of washed 
glass beads. The dashed lines correspond to the velocity data 
taken for glass beads with d = 0.5 mm and sand with d = 0.8 
mm from 



FIG. 7: (Color online) The flow velocity u as a function 
of h for a) sand with d = 0.4 mm and b) glass beads of 
d = 0.36 mm. Corresponding 6 values of 9 for a) sand 34.1° 
(o), 35.0° (a;), 36.1° (A), 37.2° (*), 40.0° (o); b) for glass 
beads 25.6° (o), 26.8° (x), 28.0° (A), 29.4° (*), 30.7° (o), 
32.0° (+). u/^/gh vs h/h s for c) sand and d) glass beads. 
u/^/gh for e) sand and f) glass beads vs. /itan 2 6/h s tan 2 6\. 



ularly for anisotropic grains we observed a larger height 
in the vicinity of the front. In other cases, typically for 
larger (spherical) grains, the front was less defined with 
some grains rolling ahead of the front, i.e., saltating. 

The collapse of the data curves for sand and glass 
beads in Figs. [71c) andJTfd) is not perfect. In these di- 
mensionless units higher plane inclinations still result in 
somewhat faster flow. We therefore consider the mod- 
ified PJFR scaling [2l[ that includes a tan 2 6 correc- 
tion to the h/h 8 term. In Figs. [TJe) andJT^f), we plot 
our data in terms of this modified scaling form, namely, 
u/y/gh versus /itan 2 9/h s tan 2 Q\. The PJFR produces 
improved scaling relative to PFR as demonstrated in 
Figs.[7Jc) andJTJd). Another consequence of the theory is 
a prediction for the density decrease with increasing 9 as 
p/p s = 1 — B-tan 6 6> (where p s stands for the close packed 
static density). Our data for the mean density, reported 
elsewhere [25[ , are well fit by the theoretical form with a 
value B = 0.52. 

We next extend our comparison of flow rules to the 
whole set of sand and glass beads used in this study. 



The data taken for these materials were scaled in the 
same manner as in the case of sand with d = 0.4 mm and 
glass beads with d = 0.36 mm as presented in Figs. [TJc) 
and[7]Jd), i.e., using PFR. In Fig. we plot uj\fgh as a 
function of h/h s for sand and glass beads with a variety 
of sizes. The scatter of the data is a sign of an imperfect 
collapse for each material. For comparison the flow rule 
measured by Pouliqucn for sand and glass beads is in- 
cluded (dashed lines) and agrees within the data scatter 
with our measurements. Plotting u/^fgh as a function 
of /itan 2 0/h s tan 2 B\ for sand and glass beads yields an 
improved collapse, see Figs. [5] and [5] 

There are two things to notice about the curves in Figs. 
|8]and[9l focusing more on the latter. First, there is the 
linearity of the lines for different d. The glass bead data 
form quite nice straight lines in support of a simple Bag- 
nold rheology with u ~ /i 3 / 2 . There is some remnant 
dependence on d discussed below. The sand data are 
quite well collapsed and have a weaker variation on the 
grain size. The curves are, however, not straight lines but 
are slightly concave downward. This deviation from lin- 
earity suggests a modification of the Bagnold rheology is 
needed but the basic form captures the main details of the 
scaling. Second, the nonzero offset 7 observed for sand 
for the case of the Pouliquen flow rule (Fig. [5]) becomes 
approximately zero for the modified scaling relationship 
(see Fig. This leads to a simpler quantitative com- 
parison of these materials as the curves are characterized 
by a single parameter (3. 
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FIG. 9: (Color online) Dimensionless flow velocity u/^fgh 
vs h tan 2 9/h B tan 2 #i for sand and glass beads as a test of 
PJFR. 



We now consider some of the details of the data from 
the perspective of the particle size d. Comparing the 
curves measured for various materials in Fig. [9l we find 
curves with similar slopes for sand of different sizes. Simi- 
larly, for glass beads the slopes (3 of the curves determined 
using PJFR and for samples of various d do not differ 
much except for the material with d = 0.72 mm where f3 
is considerably smaller. This difference cannot be quanti- 
tatively explained, but the case of the d = 0.72 mm glass 
beads could be special because the R/d ratio is very small 
(0.26) in this case. For spherical beads there is a thresh- 
old value of R/d below which the beads simply roll down 
the plane. As we approach this threshold by decreasing 
R/d the value of h s /d drops rapidly. There is a stronger 
decrease of h s /d, presumably resulting from the rolling 
effect of spherical glass beads, for spherical d = 0.72 mm 
glass beads than for irregular d = 0.85 mm sand particles 
as a function of decreasing R/d, obtained by varying the 
sandpaper roughness of R = 0.69 mm, 0.43 mm and 0.19 
mm. Using sandpaper with R = 0.12 mm the d = 0.72 
mm beads already rolled down the plane. The low value 
of h s /d could explain the low value of P measured for the 
glass beads of d — 0.72 mm on sandpaper with R = 0.19 
mm. Generally, the collapse suggests that the modified 
scaling theory describes the data quite well provided the 
roughness ratio is larger than R/d > 0.3 for glass beads 
and R/d > 0.2 for sand. 

If the flow rule provided perfect collapse of the data, 
there would be no residual dependence of the PJFR slope 
P on d. This appears to be the case for the sand flows 
where (3 w 0.37 independent of d as illustrated in Fig. 
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FIG. 10: (Color online) The PJFR slope /3 vs a) grain size d 
and b) tan 6V (x) or tan Si (o) for sand and glass bead samples. 
The dashed line corresponds to a linear fit yielding (3 = 1.22 • 
tan Si — 0.34. Slopes for copper particles (* and o) with rj = 
0.33 and 77 = 0.50 are included for comparison. 



TOTa). On the other hand, the values of (3 for glass beads 
show a systematic decrease with increasing d. It is in- 
teresting that sand with its somewhat anisotropic grains 
is less sensitive to size variation than the more idealized 
glass spheres. Again the rolling effect may play an im- 
portant role here. 

We now consider the behavior of the sand/glass-bead 
materials by plotting the slope (3 of the modified flow rule 
u/y/gh = p ■ h tan 2 9/h s tan 2 Q\ as a function of tan 9 r or 
tan 0i, see Fig. [TUT b). A significant increase in /? is ob- 
served with increasing tan 9 r or tan 0i . In a certain sense, 
these angles measure the degree of frictional interactions 
of the grains. This finding is in general agreement with 
earlier more limited data [2, H[ and gives a general char- 
acterization of the materials. Although we do not have 
enough data to unambiguously determine a functional 
dependence of P on tan0i (or tan0 r ), a linear fit to the 
data yields the relationship P = 1.22 • tan0i — 0.34. 



C. Flow rule for copper particles 

The application of the flow rule scaling to the copper 
materials is an interesting extension beyond those mate- 
rials measured previously @,[H- In particular, the copper 
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FIG. 11: (Color online) The flow velocity it as a function of h for the four set of copper particles with a) r\ — 0.25, b) rj = 0.33, 
c) rj = 0.5 and d) r\ = 0.63. e)-h) u/^/gh vs /i//i s (o) and /itan 2 #//i s tan 2 #i (x) for same rj. i)-l) ujsfgd vs /i//i s (o) and 
htan 2 9 /h a tan 2 Si (x) for same r), i)-l). 



grains are metallic and thus not affected by static charg- 
ing. Further, the grains may oxidize producing different 
frictional contacts than for the more inert sand and glass 
materials. Finally, the very unusual shape anisotropy 
adds an additional level of complexity to the scaling prob- 
lem beyond the unknown differences in shape between the 
sand and glass beads. We proceed in the same way for 
the different copper grains as with the sand/glass-bead 
materials in that first we show the raw data for u as a 
function of h in Figs. [TTT a) and [TTT d). The data vary 
smoothly with h for different values of 9. Applying PFR 
or PJFR scaling as shown in Figs. [TTT c)- fTTT h) demon- 
strates that either scaling does not work for the copper 
materials and is especially poor for the spherical copper 
grains with r\ = 0.63. The apparent origin of this poor 
collapse seems to be the assumed Y?l 2 scaling implied 
by a Bagnold vertical velocity profile. If instead of di- 
viding by y/gh, one simply plots ujsfgd versus h/h s or 
the modified form /itan 2 9/h s tan 2 9i, the curves are now 
approximately collapsed, see Figs. mTi)- fTTT l). 



In understanding this unexpected result, we first con- 
sider the spherical copper particles with r\ = 0.63 for 
which the comparison with the sand and glass beads 
might be thought to be most similar. The first thing 
to note is that there is a distinct concave downward cur- 
vature to the raw u versus h curves in Fig. Illf d) when 
compared to the case of sand or glass beads. Also, the 
character of the scaled curves is strongly non-linear for 
the case of copper with 77 = 0.63 and 77 = 0.5 (Figs. 
[TTT g) and [TTT h). Although we do not have a quantita- 
tive explanation for the behavior of the copper particle 
rheology, we note some ideas worth exploring. One issue 
of possible relevance is that the coefficient of restitution 
of soft metal particles, i.e., brass, copper, decreases with 
increasing velocity (32|, [33[ (the restitution coefficient for 
brass, which is harder than copper, decreases by about 
8% over the range of velocities in the present experiment - 
0-2.5 m/s) whereas the harder glass-bead/sand materials 
have a larger, velocity independent restitution coefficient. 
A velocity dependent (decreasing) restitution coefficient 
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would lead to higher dissipation at larger velocities, but 
this effect has not been quantitatively studied. A recent 
study on soft particles with constant restitution coeffi- 
cient [HI suggests that the presence of long-lived contacts 
leads to a modified rheology with a new term (similar to a 
Newtonian fluid), i.e., cr xz = Aj 2 + By. Such a relation- 
ship would lead to a faster growth of u with increasing ft 
than u ~ ft 3 / 2 , a result that would lead to worse agree- 
ment for our copper data than did the Bagnold scaling. 
Part of the issue here is the indirect measure of the bulk 
rheology provided by comparing the dependence of u on 
ft. 

Another possible issue is the nature of the boundary 
condition for copper particles on the sandpaper surface. 
Unlike a fluid, a granular material can have a finite slip 
velocity at the surface. This finite velocity would compli- 
cate the scaling procedure and perhaps lead to spurious 
conclusions. Copper particles move somewhat faster for 
a given thickness ft owing to their smaller size and thus 
may develop a larger slip velocity. For example, the cop- 
per particles have maximum velocities of order 2.2 m/s 
compared to 1.3-1.5 m/s for sand or glass beads over the 
same range of angle-corrected ft. 
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FIG. 12: (Color online) The dimensionless flow velocity 
u/y/gd as a function of the modified dimensionless flow thick- 
ness ft tan 2 8/h s tan 2 8i for the four set of copper particles 
with r) = 0.25, r\ = 0.33, r\ = 0.5 and 77 = 0.63. 

The other copper particles present a more complex 
situation. First, we plot the ujyfgd versus the differ- 
ent angle-corrected scalings in Fig. [T2] The data with 
?] = 0.25, 0.33, and 0.63 collapse rather well but the curve 
with r] = 0.5 has a quite different slope, about half of the 
other curves. The difference in slope does not come from 
higher velocity but rather from larger ft relative to h s . 
In other words larger h/h s was needed for the realization 
of the stationary flow regime, which results from a rel- 
atively larger dynamic friction coefficient, the source of 
which may be surface oxidation of the copper particles. 
This results in a lower value for f3. This set of copper is 
also particular in that it is the only copper sample emit- 
ting strong sound during shearing, similar to but much 
stronger than the Kelso sand. 



Although all of the copper particle data are collapsed 
better by not scaling u by yfgh, PJFR is not so bad for 
the tj — 0.33 and rj = 0.50 copper particles. Extracting a 
slope f3 for those values of 77 yields curves that are con- 
sistent with the sand/glass-bead scaling as a function of 
tan 6*i, see Fig. [TO] Thus, even though the copper parti- 
cles are quite different, they still seem to show the same 
qualitative dependence on tan#i as the glass-bead/sand 
particles. 



IV. CONCLUSIONS 

The most important findings of this work can be sum- 
marized as follows. The surface velocity u as a func- 
tion of flow thickness ft, of a granular flow on a rough 
inclined plane was measured for fourteen different ma- 
terials in the dense, stationary flow regime. All config- 
urations were characterized by measuring the value of 
h s (the thickness of the layer remaining on the plane 
after the flow subsided) as a function of the plane incli- 
nation 9. The value of h s /d for sand and glass beads 
increased slightly with increasing ratio of plane rough- 
ness and grain diameter R/d measured for four different 
values of R. The u(h) curves for sand and glass beads 
measured at various 9 did not perfectly collapse using 
the scaling law u/y/gh ~ h/h s proposed by Pouliquen 
0. An improved collapse was obtained using PJFR 
u/y/gK = P ' ht&n 2 9/h s tan 2 9\ where the factor tan 2 9 
was suggested by a recent theory by Jenkins [2l[ . For the 
sand/glass-bead materials, the PJFR slope /3 increases 
strongly with tan 9\ yielding a quantitative description of 
various materials, thereby extending our tools for a bet- 
ter characterization and prediction of complex dynamical 
phenomena, such as waves 0] or avalanche propagation 

Our results demonstrate that when the surface velocity 
is used to determine the flow rule, the PJFR scaling is su- 
perior to the earlier PFR approach. For the original data 
set of Pouliquen , it is hard to determine which scaling 
form is better. Two possibilities are suggested. First, the 
uncertainty (mainly in h s ) in the original measurements 
does not allow a definitive comparison. Second, measur- 
ing the depth averaged velocity at the front is substan- 
tially different from measuring the surface velocity. If the 
latter is correct then something unexpected is happening 
in the layer because Bagnold scaling (or any monotonic 
vertical velocity profile starting from zero velocity, for 
that matter) implies that the ratio of surface velocity to 
depth averaged velocity is a constant, and a constant fac- 
tor would not change the flow rule. Although we cannot 
definitively rule out some strange behavior, the superior 
fit of PJFR, the elimination of the need for an offset 7, 
and the consistency of these results with Bagnold scal- 
ing suggests that the apparent discrepancy between our 
results and former scaling analysis 0] results from larger 
uncertainty in the previous measurements. 

For copper grains of different shapes, neither the 
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Pouliqucn form nor the Jenkins modified scaling works 
well in collapsing data taken for a variety of values of 9. 
Although the angle correction works decently, the nor- 
malization of u by \fgh produces poorer scaling. This 
suggests that for the copper particles, a Bagnold form 
for the vertical velocity profile does not hold. An impor- 
tant future extension of this work would be to directly 
measure the velocity as a function of vertical position for 
the different materials to determine the velocity profiles. 
Measurements of this type are being planned to test the 
conjectures based on phenomenological flow rule compar- 
isons. 

Finally, one must conclude that although the rheology 
for sand and glass beads seems rather robust and well 
fit by the Pouliqucn/ Jenkins form, this is no guarantee 
that more general materials satisfy this scaling relation- 
ship. The copper measurements are puzzling because one 



might have expected the nearly spherical copper beads 
to produce results similar to the spherical glass beads. 
That the Bagnold form does not seem to apply for copper 
grains of different shapes and especially for the spherical 
ones is quite surprising and unexpected. Experiments on 
other metallic particles would be very helpful in deter- 
mining the origins of this effect. Finally, a more direct 
probe of the interior dynamics of granular flows seems es- 
sential for determining the bulk flow rheology for general 
granular media. 
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